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ABSTRACT 


INTRODUCTION 


1.  Scope  of  the  thesis 

— -1  Careful  experiments  on  boundary  layer  transition  conducted  in  low  tur¬ 
bulence  wind  tunnels  have  resolved  three  stages  in  the  breakdown  of  lam¬ 
inar  flow  into  turbulence:  in  the  first  stage,  two-dimensional  oscillations, 
commonly  referred  to  as  TS  waves,  appear  and  propagate  in  the  direction  of 
the  flow;  in  the  second  stage  there  appears  a  rapid  three-dimensional  defor¬ 
mation  of  the  oscillatory  flow;  and  in  the  brief  third  stage  high  frequency, 
small  scale,  quasi-random  fluctuations  appear  and  quickly  lead  the  flow  into 
turbulence. 

The  recent  analysis  of  the  second  stage  based  on  a  secondary  instability 
mechanism  arising  from  a  parametric  resonance  induced  by  the  TS  wave  has 
given  results  in  quantitative  agreement  with  experimental  observations. 
Using  this  approach,  Herbert  has  studied  the  phenomenon  in  plane  channel 

h  '  '  •  - 

flow  and  in  a  Blasius  boundary  layer  [1-5].  The  present-work  studies  the 
effect  of  pressure  gradients  on  the  behavior  of  secondary  instabilities,  by 
extending  Herbert’s  analysis  to  the  Falkner-Skan  family  of  profiles,  and 
develops  a  spatial  formulation  for  disturbance  growth  to  exactly  match  the 
experimentally  observed  behavior.  L  rp,< ? 


2.  Motivation 


The  study  of  secondary  instabilities  in  the  Falkner-Skan  family  of 
profiles,  and  hence  in  flows  with  pressure  gradients,  is  a  necessary  step 
towards  the  analysis  of  transition  on  objects  with  a  non-constant  surface 
pressure  distribution,  foremost  among  which  are  airfoils.  In  high- 
Reynolds-number  aerodynamics  the  effect  of  transition  is  mainly  a  quantita¬ 
tive  change  in  drag.  At  low  Reynolds  number,  and  in  combination  with  a 
separating  boundary  layer  in  strong  adverse  pressure  gradients,  the  location 
of  transition  and  the  structure  of  the  transitional  flow  can  affect  the  global 
flow  field  about  the  airfoil.  Indeed,  with  the  presence  of  the  second  stage  of 
transition  in  a  boundary  layer  undergoing  separation,  the  three-dimensional 
disturbances  can  activate  strong  vortex  stretching  in  the  subsequent  free 
shear  layer  and  cause  breakdown  into  wall  bound  turbulence,  forming  a 
separation  bubble,  while  the  absence  of  three  dimensionality  at  separation 
can  allow  periodical  vortex  pairing  in  the  free  shear  layer  and  postpone  tur¬ 
bulence  far  enough  downstream  to  cause  the  layer  to  permanently  leave  the 
surface,  creating  a  condition  of  stall  over  the  airfoil.  A  visual  representation 
of  this  idea  is  given  in  figure  1.  Hence  in  both  high  and  low  Reynolds 
number  external  flows  the  presence  and  structure  of  transition  is  an  impor¬ 
tant  phenomenon. 

In  the  past  eight  decades,  starting  with  the  work  of  Orr  and  Sommer- 
feld,  extensive  study  has  been  done  on  the  first  stage  (primary  instabilities) 
using  the  method  of  small  oscillations.  This  approach  yields  correctly  the 
downstream  location  (critical  Reynolds  number)  at  which  instability  ini¬ 
tiates,  and  describes  accurately  the  oscillatory  field  at  small  amplitudes. 


However,  it  cannot  model  the  evolution  of  the  three-dimensional  stage, 
and,  when  used  in  conjunction  with  a  en  transition  criterion,  it  gives  transi¬ 
tional  predictions  erroneously  independent  of  the  upstream  background  dis¬ 
turbance  amplitudes. 

The  secondary  instability  model  used  in  the  present  work  overcomes 
these  shortcomings  by  incorporating  the  background  disturbance  amplitudes 
as  a  parameter,  and  predicting  in  agreement  with  observation  the  form  and 
growth  rate  of  secondary  disturbances.  We  should  mention,  however,  that 
this  model  is  not  complete  in  itself,  for  it  does  not  include  a  criterion  for 
determining  the  conditions  at  which  the  three-dimensional  disturbances 
develop  a  “life  of  their  own”  by  becoming  independent  from  the  parametric 
excitation  of  the  TS  wave  as  a  source  of  growth.  This  step,  along  with  the 
usage  of  actual  airfoil  boundary  layer  profiles  in  a  downstream  marching 
transition  prediction  analysis,  will  be  a  future  extension  of  the  present  work. 

3.  Experimental  view  of  boundary  layer  transition 

In  the  study  of  boundary  layer  transition  the  majority  of  the  experimen¬ 
tal  results  have  been  obtained  by  using  the  vibrating  ribbon  technique  of 
Schubauer  &  Skramstad  [6],  in  which  the  transition  over  a  body  immersed 
in  a  low  turbulence  stream  is  initiated,  and  hence,  controlled  by  the  regu¬ 
lated  vibration  of  a  tightly  stretched  ribbon.  A  more  detailed  account  of 
the  following  descriptions  can  be  found  in  Klebanoff,  Tidstrom,  Sc  Sargent 
[7],  Kachanov  Sc  Levchenko  [8]  ,  or  in  Saric  et  al.  [9].  Let  the  streamwise, 
normal-to-plate,  and  spanwise  directions  compose  an  orthogonal  system.  In 
a  low  noise  environment  a  streakline  observation  in  a  boundary  layer  con- 
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taining  a  forced  oscillation,  introduced  by  a  vibrating  ribbon,  might  proceed 
as  follows: 

•  At  low  oscillation  amplitudes  one  observes  traveling  waves  (TS  waves) 
invariant  in  the  spanwise  direction  and  propagating  with  the  the  flow. 
Their  amplitude  harmlessly  grows  and  decays  in  full  accordance  with 
linear  stability  theory.  Turbulence  is  delayed  to  high  Renoylds 
numbers. 

•  When  the  amplitude  of  the  TS  wave  is  increased  and  a  certain  critical 
combination  of  Reynolds  number  and  amplitude  is  reached,  the  wave 
experiences  a  rapid  spanwise  periodic  deformation.  One  can  observe 
spanwise  alternating  peaks  and  valleys  which  indicate  the  presence  of 
enhanced  and  reduced  wave  amplitude.  The  peak-valley  structure  is 
staggered,  repeating  itself  with  twice  the  periodicity  of  the  TS  wave  (see 
Fig.  2a),  and  fixed  hot  wire  measurements  record  subharmonic  signals. 
The  rapid  growth  of  the  three-dimensional  structure  leads  to  the  formar 
tion  of  concentrated  shear  layers  at  the  peak  positions.  The  highly 
inflectional  velocity  distribution  quickly  breaks  up  into  tertiary  instabili¬ 
ties. 

•  At  higher  initial  levels  of  the  TS  amplitude  the  above  phenomenon  is 
repeated,  but  with  a  different  structure;  the  peaks  and  valleys  are 
aligned,  repeating  themselves  with  the  same  frequency  and  periodicity 
as  the  TS  wave.  A  system  of  streamwise  vortices  appears  simultane¬ 
ously  with  the  aligned  peaks  and  valleys  (see  Fig.  2b).  Being  first 
observed  in  1962  by  Klebanoff  et  al.  it  is  usual  today  to  name  this 


aligned  peak-valley  splitting  as  K-type  breakdown.  The  disturbances 
undergo  pure  spatial  growth  in  all  the  above  cases. 

4.  Review  of  past  research 

The  three  stages  of  boundary  layer  transition  described  have  also  been 
observed  in  plane  Poiseuille  flow  [10],  in  free  shear  layers  [11],  and  in  a 
fluid  with  continuously  stratified  density  (internal  gravity  waves)  [12].  In  all 
three  flows,  researchers  have  tried  to  model  the  second  stage  of  transition 
using  two  approaches:  weakly  nonlinear  models,  or  parametric  instability 
models.  The  latter  one  has  given  the  most  general  and  accurate  results. 

Weakly  nonlinear  models  appeared  following  the  innovative  work  of 
Landau.  As  discussed  by  Drazin  and  Reid  [13],  one  can  expand  the  velocity 
field  in  weakly  nonlinear  analysis  as 

v  =  A,-(  t)ft-(x)  +  Complex  Conj.  (1) 

»  =l 

where  {  f,  }  is  a  complete  set  of  functions  satisfying  the  boundary  condi¬ 
tions,  and  is  usually  taken  to  be  the  wave  type  modes  of  the  linear  problem 
(Orr-Sommerfeld  or  Squire  modes  in  boundary  layers,  Helmholtz  modes  in 
free  shear  layers).  When  the  series  is  inserted  into  the  Navier-Stokes  equar 
tion  the  nonlinear  convective  term  (vv)v  creates  resonant  wave  interac¬ 
tions.  The  interaction  is,  by  computational  necessity,  limited  to  a  finite 
number  of  waves  by  truncating  the  series  to  a  few  terms,  usually  three. 
Different  weakly  nonlinear  models  stem  from  different  selections  of  modes 
in  the  truncated  series.  For  compactness  we  introduce  the  wave  vector 
k  =  (a,/?),  where  a  is  the  streamwise  wavenumber  and  3  is  the 


(perpendicular)  spanwise  wave  number. 

The  model  of  Benney  and  Lin  [14],  1960,  attempts  to  describe  the 
Klebanoff  type  structure  by  considering  the  interaction  between  the  TS 
wave,  k0  =  (a,0),  and  two  oblique  modes  of  the  Orr-Sommerfeld  equation 
(OSE),  kj  =  (a,  +0)  and  kj  —  (a,  -0),  with  0  as  a  free  parameter. 
Whether  or  not  conditions  for  phase  synchronization  of  the  three  modes 
exist  is  still  a  subject  of  debate  [15].  The  model  due  to  Craik  [16]  uses  the 
TS  wave  and  two  oblique  modes  of  the  OSE,  kj  =  (a/2,  +0*)  and 
k2  =  (a/ 2,  -  0*)  to  model  the  subharmonic  structure.  The  value  of  0*  is 
chosen  such  that  the  modes  ki  and  k2  have  the  same  phase  speed  as  the  TS 
wave,  thus  allowing  for  efficient  energy  transfer  among  the  three  modes.  In 
boundary  layers,  this  model  sometimes  predicts  correctly  the  value  spanwise 
undulation  0,  but  at  other  times,  under  similar  experimental  conditions,  0* 
differs  widely  from  the  observed  value.  In  plane  Poiseuille  flow  the  model 
is  inoperative  for  reasons  of  symmetry  [17].  Herbert  and  Morkovin  [17]  for¬ 
warded  a  model  composed  of  the  TS  wave  and  two  longitudinal  vortices  of 
the  OSE  equation,  ^  =  (0,  +0)  and  k2  =  (0,  -  0),  to  model  a  spanwise 
periodic  longitudinal  vortex  structure. 

Other  models  have  been  developed,  and  the  above  list  is  by  no  means 
exhaustive.  However,  along  with  the  weakly  nonlinear  model  the  author 
should  present  a  guarantee  of  “completeness”  of  the  model,  i.e.  show  that 
the  modes  considered  are  responsible  for  the  observed  phenomena.  This 
task  is  difficult,  it  is  rarely  performed,  and  is  the  single  greatest  drawback  of 
the  method. 

The  secondary  instability  model  based  on  parametric  resonance  is 
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intuitively  simple  to  understand  using  physical  arguments.  In  the  first  stage 
of  transition  the  presence  of  the  primary  wave  brings  about  a  redistribution 
of  vorticity  in  the  flow.  When  moving  with  the  velocity  of  the  wave,  an 
observer  sees  a  periodically  spaced,  coplanar  arrangement  of  vortex  tubes, 
sometimes  referred  to  as  a  “cat’s  eyes”  vorticity  pattern.  This  arrangement 
is  highly  unstable  to  small  three-dimensional  disturbances:  a  vortex  tube 
slightly  distorted  along  its  length  is  alternatively  placed  into  regions  of  faster 
and  slower  moving  fluid  and  undergoes  rapid  stretching.  Furthermore,  if 
there  is  sufficient  room,  the  arrangement  is  also  unstable  to  small  two- 
dimensional  disturbances  which  tend  to  amalgamate  neighboring  vortices 
into  a  larger  vortex  (vortex  pairing).  The  mathematical  analysis,  hence,  can 
be  divided  into  two  parts;  the  construction  of  the  the  periodic  flow  arising 
from  the  first  stage  of  transition,  and  the  analysis  cf  a  superimposed  three- 
dimensional  velocity  perturbation  of  amplitude  sufficiently  small  for  lineari¬ 
zation.  Due  to  the  spatial  periodicity  of  the  basic  flow,  obe  instability 
mechanism  is  governed  by  a  Floquet  system. 

The  first  application  of  a  Floquet  system  to  boundary  layer  transition 
seems  to  have  been  done  by  Maseev  [18],  1968.  However,  his  work  went 
unobserved,  perhaps  due  to  its  short  presentation  and  use  of  an  unclear  per¬ 
turbation  method  for  solution.  Prior  to  this  work,  Floquet  theory  had  been 
used  by  Kelly  [19],  1967,  for  the  stability  study  of  a  inviscid  free  shear 
layer.  He  found  the  vortex  pairing  mode  to  stem  from  principal  parametric 
(subharmonic)  resonance.  Pierrehumbert  &  Wldnall  [20].  1979.  repeated 
the  analysis  using  a  more  accurate  periodic  flow  model  composed  of  an  array 
of  Stuart  vortices. 


Observing  that  certain  experimental  setups  for  the  study  of  Lamb- 
Taylor  instabilities  excited  periodically  the  flow,  Nayfeh  [21],  1969,  derived 
a  Mathieu  type  equation  to  explain  the  overstability  behavior.  In  the  mid 
and  late  seventies  several  researchers  applied  Floquet  theory  to  model  the 
evolution  of  three-dimensionality  in  internal  gravity  waves.  McEwan  & 
Robinson  [12],  1975,  studied  parametric  instability  driven  by  time  periodic 
coefficients.  Mied  [22],  1976,  refined  the  analysis  by  changing  to  a  frame 
moving  with  the  periodic  wave  and  solved  for  temporally  growing  distur¬ 
bances,  and  extended  the  analysis  to  Rossby  waves  [23].  In  both  studies  he 
shows  that  “the  parametric  instabilities  reduce  to  the  nonlinear  resonant 
interaction  in  the  limit  of  vanishing  basic  state  amplitude’’.  Drazin  [24], 
1977,  independently  developed  a  Floquet  formulation  for  this  problem,  and 
in  his  work  he  related  the  parametric  instability  to  resonant  wave  interaction 
and  to  the  catastrophe  theory  of  Thom.  He  noted  that  the  Floquet  system 
could  be  applied  to  other  flows  if  certain  general  properties  of  waves  are 
satisfied. 

In  1979,  Nayfeh  and  Bozatli  [25]  studied  the  subharmonic  vortex  pair¬ 
ing  in  a  boundary  layer  using  the  method  of  multiple  scales.  They  showed 
that  a  large  amplitude  of  the  TS  wave  (29%)  is  necessary  to  initiate  the  pair¬ 
ing  process.  This  fact  is  in  agreement  with  the  present  results.  However  the 
analysis  was  only  two-dimensional,  while  the  most  unstable  subharmonic 
modes  are  three-dimensional,  with  a  spanwise  wave  number  of  the  order  of 
the  TS  wave. 

The  application  of  Floquet  theory  to  the  stability  study  of  equilibrium 
states  in  plane  Poiseuille  flow  was  concurrently  developed  by  Orszag  [26] 


and  Herbert  [27],  1981.  The  former  author  used  the  method  to  verify  the 
linear  character  of  three-dimensional  instability  as  revealed  by  the  distur¬ 
bances’  exponential  growth  with  time  in  a  full  numerical  simulation  of  the 
Navier-Stokes  equation.  While  Orszag  focused  on  large  amplitude  equili¬ 
brium  states  and  subcritical  instability,  Herbert  performed  a  thorough  inves¬ 
tigation  in  the  case  of  small  amplitude  TS  waves,  and  used  the  threshold 
amplitudes,  the  /?  selectivity,  and  other  characteristics  of  secondary  instabili¬ 
ties  to  clarify  the  experimental  observations.  He  extended  the  model  to  the 
Blasius  boundary  layer.  By  a  careful  study  in  the  limit  of  vanishing  TS  wave 
amplitude  he  was  able  to  determine  some  circumstances  in  which  weakly 
nonlinear  models  were  complete  and  others  in  which  the  models  failed  to 
account  for  the  dominant  modes  [4].  By  means  of  a  crude  transformation 
from  temporal  to  spatial  growth  rates,  the  experimentally  observed  stream- 
wise  amplitude  variation  of  the  subharmonic  was  astonishingly  reproduced 
by  the  theory  [2],  confirming  the  parametric  nature  of  instability  in  the 
second  stage  of  transition. 


GOVERNING  EQUATIONS  FOR  SECONDARY 


INSTABILITY 


The  following  derivation  is  confined  to  steady,  incompressible,  constant 
temperature  shear-layer  flows  having  a  weak  variation  in  the  streamwise 
coordinate.  In  the  derivation  we  follow  Herbert  [1-5].  The  analysis  is 
divided  into  two  parts:  the  construction  of  the  periodic  basic  flow  and  the 
investigation  of  superposed  secondary  disturbances  with  amplitudes 
sufficiently  small  for  linearization. 


1.  Construction  of  the  Periodic  Basic  Flow 

We  consider  the  flow  over  a  semi-infinite  plate  at  an  angle  P tt/2  to  the 
free  stream.  Let  X  denote  the  distance  from  the  leading  edge  along  the 
plate,  Y  the  distance  normal  to  the  plate,  and  z  the  spanwise  coordinate, 
U,  =  kXp^2~p'>  the  inviscid  flow  velocity  at  the  surface,  and  X0  the  loca¬ 
tion  at  which  the  stability  analysis  is  performed.  The  undisturbed  boundary 
layer  stream  function  is  given  by 


*FS(X,Y,t)  =  U,f(,)S  ,  S  =  ,  y 

where  / (y)  satisfies  the  Falkner-Skan  equation 

/"'  +  U"  +  P{  l-  f2)  =  0 
/( 0)  -  /'CO)  =  0  ,  /'(y-co)  =  1 
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We  use  <5  and  Ut  as  characteristic  quantities  for  nondimensionalization,  and 
define  the  Reynolds  number  by  R  =  U,6/ u.  At  sufficiently  high  R,  pri¬ 
mary  disturbances  appear,  and  during  their  initial  stages  of  growth  they  are 
well  modeled  by  linear  stability  analysis.  Only  a  general  overview  of  the 
analysis  will  be  given  here.  For  a  detailed  description  see  ,for  example, 
Betchov  and  Criminale  [28].  We  superimpose  a  small  disturbance  vx  on  the 
basic  flow  V, 

v(X,y,z,t)  =  V(X,y)  +  A  vv(X,y,z,t)  (4) 

A  «  1  (5) 

and  insert  this  compound  flow  into  the  Navier- Stokes  equation.  The  non¬ 

linear  convective  term  A2(v1v)  is  dropped  in  view  of  the  small  ampli¬ 
tude  of  the  disturbance.  Collecting  terms  of  order  A  one  obtains  linear 
governing  equations  for  vx, 

■j£y>  +  (vrv)V  +  (V-v)vi  =  -  vj>i  +  jt-v2v,  (6) 

W,  =  o  (7) 

subject  to  the  boundary  conditions  v^XjO,.?,*)  =  0  , 

vr(X,  y— *oo  ,z,t )  =  0.  Equation  (6)  is  separable  in  t  and  z.  Due  to  the 
weaJc  variation  of  the  Falkner-Skan  profile  with  X,  we  apply  the  parallel-flow 
approximation  and  assume  the  profile  to  be  locally  one-dimensional, 
V0(X9,y).  Then  the  disturbance  equation  becomes  separable  also  in  X,  and 
the  solution  can  be  expressed  with  normal  modes  in  the  variables  X  ,  t  , 
and  z.  The  error  associated  with  the  parallel-flow  approximation  is  small  in 
profiles  with  favorable  pressure  gradients,  P  >  0,  but  increases  in  profiles 
with  adverse  pressure  gradients,  P  <  0.  For  a  detailed  description  of  the 


effect  of  non-parallelism  in  the  analysis,  see  Saric  &  Nayfeh  [29].  As  first 
shown  by  Squires  [30],  the  critical  Reynolds  number  is  associated  with  a 
two-dimensional  wave.  This  fact  has  lead  to  the  incorrect  popular  belief  that 
two  dimensional  disturbances  are  most  unstable  at  all  values  of  wavenumber 
a,  or  frequency  F.  In  our  analysis,  we  restrict  our  attention  to  two- 
dimensional  TS  waves  since  this  type  is  most  easily  produced  by  the  vibrat¬ 
ing  ribbon  technique.  The  usage  of  oblique  TS  waves  in  our  analysis  should 
require  little  modification.  Thus,  we  drop  the  z  dependence  and  recast  the 
equations  in  terms  of  the  disturbance  stream  function  The  solution  form 
is  given  by 

=  A( X)  [*(?)«'“'<*- “>  +  *>(,),'  (8) 

The  t  denotes  the  complex  conjugate,  c  is  the  phase  velocity,  a  =  a:r+ta,- 
and  <t>(y)  are  solutions  to  the  Orr-Sommerfeld  eigenvalue  problem, 

{D2-iaR((U,  -  c)D  -  U /))*  =  0  (9) 

0(0)  =  <p'(0)  =  0  ,  0(-*oo)  =  0'(-*oo)  =  0 

d2 

where  D  =  — r  -  ex2  and  ar  is  the  streamwise  wavenumber,  a{  is  the 
dy 2 

growth  rate.  Note  that  R  is  a  coefficient  in  (9),  so  that  a  =  a  (R)  and, 
consequently,  a  is  also  a  function  of  the  streamwise  position  X.  The  variar 
tion  of  the  amplitude  is  given  by 

X 

A(X)  =  A,,exp(-/  a,(X)dX)  (10) 

■V. 

where  A,„  is  the  TS  amplitude  at  an  initial  location  In  the  presence  of  a 
TS  wave  of  finite  amplitude,  the  two-dimensional  flow  has  the  form 


v2(X,y,t)  =  V0(y)  +  A(X)  v^X-  ct,y ) 


(11) 


The  eigenfunction  <j>{y)  in  (8)  is  normalized  such  that  A  measures  the 
streamwise  rms  fluctuation  of  the  disturbance.  As  the  primary  disturbance 
grows  from  infinitesimal  to  finite  amplitude,  nonlinear  interactions  begin  to 
modify  both  the  mean  and  the  disturbance  profiles.  However,  the  numerical 
studies  of  Orszag  [26]  and  Herbert  [4]  in  plane  Poiseuille  flow  have  shown 
that  nonlinear  modifications  have  negligible  effect  on  the  growth  of  secon¬ 
dary  instabilities.  Thus,  we  keep  the  shape  of  the  TS  wave  as  given  by  the 
linear  analysis  up  through  finite  amplitude  values.  This  approximation  is 
called  the  shape  assumption. 

In  view  of  the  weak  (viscous)  growth  of  the  TS  wave  in  comparison  to 
the  explosive  (convective)  growth  of  secondary  disturbances,  we  neglect  the 
weak  variation  of  the  TS  amplitude  and  consider,  locally, 
A(X)  —  A(X0)  =  A0.  One  may  view  this  approximation  as  an  extension  of 
the  parallel  flow  assumption  to  the  TS  wave,  and  as  the  order  approxi¬ 
mation  in  a  weakly  non-parallel  analysis.  The  time  dependence  is  eliminated 
by  changing  from  a  fixed  coordinate  system  to  a  coordinate  system  moving 
with  the  phase  velocity  of  the  TS  wave,  x  =  X  -  ct  ,  and  the  flow 
becomes  periodic  in  the  moving  streamwise  coordinate  x, 

v2(i,y)  =  V9(y)  +  A9vt(x,y)  (12) 

Here,  and  elsewhere,  the  'symbolizes  the  periodicity  of  the  function. 

Fig.  3  illustrates  the  above  steps.  In  a  neighborhood  of  X0  both  the 
Falkner-Skan  profile  and  the  TS  wave  amplitude  are  assumed  constant, 
thereby  extending  the  periodic  flow  from  -  cc  to  oc.  A  detailed  study  of 


-  14  - 


temporally  and  spatially  growing  TS  waves  for  Falkner-Skan  profiles  has 
been  presented  by  Wazzan  et  al.  [31].  The  periodic  basic  flow  (12)  can  be 
considered  as  known  from  any  standard  procedure  for  analysis  of  TS  waves. 


2.  Secondary  Disturbances.  Linear  model. 

We  study  the  evolution  of  secondary  disturbances  through  a  linear  star 
bility  analysis  of  the  periodic  basic  flow  (12).  A  small  secondary  disturbance 
v3  —  (u3  >v3  >w3,)  is  superimposed  on  the  basic  flow, 

y{x,y,z,t)  —  v2(x,y)  +  ev3(x,y,z,t)  ,  e  «  1  (13) 

and  the  compound  flow  is  inserted  into  the  Navier-Stokes  equations.  The 
nonlinear  convective  term  e2  (v3-^)v3  is  dropped  in  view  of  the  small  dis¬ 
turbance  amplitude,  yielding  linear  equations  for  v3, 

d  1 

-£jv3  +  (v2  v)v3  +  (v3-v)v2  *  -VPs  +  ^*V2v3  (14a) 

vv3  =  0  (14b) 


subject  to  the  boundary  conditions  v3(x,0,2,£)  =  0  , 

v3(x,  y— *-oo  ,z,t)  =s  0.  The  linear  secondary  stability  problem,  hence,  is  an 
eigenvalue  problem.  Taking  the  curl  of  (14a)  to  eliminate  the  pressure, 
introducing  =  ~V2^i  and  the  disturbance  vorticity  u/3  =  (£3  ,r?3  ,^3), 
and  using  continuity  to  eliminate  the  >component  of  velocity  tv3,  we  obtain 
a  system  of  of  linear  differential  equations  for  u3  and  v3, 


1  2 
"TTv  - 


f  rr  ,3  9 

"  C)  dx  '  dt 


V*v3 


+ 


d2u,  dv3 

dy 2  dx 


ggi  d 

dx  dy 


drj>  1  d 
dy  dx 
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Note  that  dt]3fdz  and  df3 jdz  are  expressible  in  terms  of  u3  and  u3  .  These 
equations  have  three  important  features:  linearity,  coefficients  independent 
of  time  and  spanwise  coordinate  z,  and  coefficients  periodic  in  x  due  to  the 
presence  of  vt.  The  first  two  allow  the  use  of  normal  modes  in  t  and  z, 

v3(x,y,*,f)  =  e,0ze<Tt  v3(x ,y)  (16) 

where  the  spanwise  wavenumber  0  is  real,  and  a  is  in  general  complex. 
Floquet  theory  is  fully  developed  in  the  context  of  ordinary  differential 
equations,  and  is  well  documented  in  the  literature  (see  for  example  Nayfeh 
&  Mook  [32],  or  Arscott  [33]).  The  extension  to  partial  differential  equa¬ 
tions  requires  sophisticated  mathematics,  but  seems  to  provide  very  similar 
basic  results.  In  our  problem  we  assume  a  solution  of  Floquet  form,  since, 
apart  from  the  y  dependence,  the  secondary  disturbance  equations  are 
essentially  of  Hill  type  in  x.  Letting  7  denote  the  characteristic  exponent, 
and  X,  =2 tt /  ar  the  TS  wavelength,  we  find  two  classes  of  solutions  v3  in 
the  form 

v3(x,y)  =  e7*  v/(*,y)  v/(x  +  Xx,y)  =  v/(x,y)  (17a) 

v3(x,y)  =  e71  y3(x,y)  v3(x  +  2Xz,y)  =  vt(x,y)  (17b) 

where  we  denote  Vy(x,y)  as  the  “fundamental”  mode,  and  v3(x,y)  as  the 
“subharmonic”  mode.  Being  periodic,  Vy  and  v,  can  be  represented  by 
Fourier  series. 


*f(x,y)  =  £  v2n(y)e‘2nd*  ^  (18a) 

»  =-  oo  a  =  a/2 

♦.(*.»)-  E  ^.  +  1(Sf)«'(2"+l)d'  (18b) 

n  =  -  oo 

In  principle,  we  can  arrive  at  the  same  result  by  first  expressing  the  solution 
in  a  Fourier  series  in  a, 

g  v.tsiie”6'  (19) 

n  =*-  oo 

and  then  introducing  the  basic  flow  (12)  and  secondary  disturbance  (19) 
into  the  equations  (15).  The  terms  multiplied  by  exponentials  of  like  index 
are  collected  and  set  equal  to  zero,  two  uncoupled  infinite  sets  of  equations 
result  for  the  fundamental  components  v2n(y)  and  subharmonic  com¬ 
ponents  v2n  +  1(y).  Due  to  the  extensiveness  of  algebra  involved,  these 
equations  have  been  relegated  to  the  appendix. 

The  solution  forms  (17a,b)  reproduce  the  experimentally  observed 
behavior  of  three-dimensional  disturbances  [1,3].  The  fundamental  mode 
produces  the  Xt-periodic  “in-line”  or  peak-valley  splitting  pattern  as 
observed  by  Saric  &  Thomas  [9,  Fig.  l].  The  aperiodic  term  v0(y)  in  the 
series  (17a)  gives  rise  to  both  a  mean  flow  distortion  and  a  spanwise 
periodic  longitudinal  vortex  structure.  The  subharmonic  mode  produces  the 
2Xj-periodic  staggered  structure  [9,  Fig.  2,  3]  and  has  no  aperiodic  term  in 
the  series  (17b).  Both  the  fundamental  and  subharmonic  modes  can 
undergo  exponential  growth  in  time  or  space,  consistent  with  experience 
from  numerical  simulations  and  experiments  (see  also  [34]). 

The  secondary  instability  equations  can  be  recasted  with  real 
coefficients,  since  the  Navier-Stokes  equation  involves  only  real  values. 


Thus,  the  eigenvalues  are  either  real,  or  complex  conjugates.  In  the  case  of 
complex  eigenvalues  special  care  must  be  taken  in  the  computational  tech¬ 
nique  to  follow  the  birth  and  collapse  of  the  complex  parts.  For  real  eigen¬ 
values,  the  method  of  solution  herein  presented  is  valid,  and  the  use  of  a 
real  formulation  provides  considerable  savings  in  the  numerical  work. 

3.  Spatial  and  temporal  growth  in  a  moving  frame 

Owing  to  the  exponential  dependence  of  Vj  on  both  time  and  space,  as 
seen  in  (19),  we  are  faced  with  the  necessity  of  choosing  one  of  the 
exponents  ,  a  or  7,  as  the  eigenvalue.  In  this  manner  we  select  between 
temporally  growing  or  spatially  growing  modes,  in  full  similarity  to  the  study 
of  TS  waves.  Temporal  growth  analysis  is  simpler  and  requires  less  compu¬ 
tational  effort  since  a  appears  linearly  in  the  equations,  while  7  appears  up 
to  fourth  power.  Moreover,  a  relation  exists  between  the  two  types  of 
growth,  and  can  be  utilized  to  approximately  transform  temporal  into  spatial 
growth  rates  (see  below). 

For  primary  disturbances,  spatial  growth  is  measured  with  respect  to  the 
laboratory  frame  of  coordinates.  In  order  to  arrive  at  a  consistent  formula¬ 
tion  of  temporal  and  spatial  amplification  for  secondary  modes,  we  express 
the  velocity  field  (19)  in  laboratory  coordinates 

v3(X,y,z,t)  =  £  vn{y)t,n6(X-  ct)  (20) 

n  =-  00 

To  restrict  the  disturbance  exclusively  to  temporal  amplification  we  set 
the  “spatial”  exponent  7  to  zero  and  solve  for  the  “temporal”  exponent  a 
as  an  eigenvalue.  The  real  part  of  a  dictates  growth  and  decay  in  both 


coordinate  systems,  hence  the  moving  and  the  laboratory- fixed  observer  wit¬ 
ness  the  same  temporal  growth.  For  purely  real  a ,  the  mode  is  fixed  in  the 
( x,y,z,t )  frame  and,  hence,  tuned  with  the  TS  wave.  For  complex  a  one 
obtains  conjugate  pairs,  a  =  o  T±  a and  the  modes  propagates  upstream 
and  downstream  in  the  moving  coordinate  system.  In  a  laboratory  fixed 
frame  the  disturbances  travel  with  a  phase  speed  different  from  that  of  the 
TS  wave.  Calculations  have  shown  that  at  realistic  amplitudes  the  dom¬ 
inant,  i.  e.  most  unstable,  mode  is  associated  with  a  real  eigenvalue. 

A  spatially  amplified  disturbance  in  the  moving  coordinate  system  can 
be  computed  by  setting  a  equal  to  zero  and  solving  for  7  as  an  eigenvalue. 
To  a  lab  oratory- fixed  observer,  however,  this  disturbance  mode  exhibits 
mixed  spatial-temporal  amplification,  since  in  (15)  both  a  and  7  are  multi¬ 
plied  by  time. 

The  two  cases  of  most  interest  are  those  of  either  temporal  or  spatial 
amplification  in  laboratory  coordinates.  Indeed,  experimental  observation 
shows  that  in  boundary  layers  the  growth  of  primary  and  secondary  distur¬ 
bances  is  only  spatial. 

In  order  to  obtain  disturbance  modes  with  purely  spatial  amplification  in 
the  laboratory  frame,  we  note  from  (20)  that  the  time-dependent  factor 
e((T-'7c)t  can  be  suppressed  by  choosing  a  =70.  In  the  moving  frame,  we 
therefore  search  for  solutions  of  the  form 

v3(x,y,z,0  =  e^e^e^  £  vn(y)e"*«  (21) 

n  =1-  00 

and  solve  for  7  as  the  eigenvalue.  (Recall  that  computations  must  be  carried 
out  in  the  moving  coordinate  system  for  only  there  is  the  Floquet  theory 


directly  applicable).  We  calculated  all  the  information  on  spatially  growing 
modes  in  this  report  using  this  approach.  Lastly,  we  note  that  setting 
a  =  ~tc  +  a  produces  disturbances  growing  spatially  in  X  but  having  the 
frequency  detuned  from  the  subharmonic  by  the  value  of  d.  In  this  case  the 
eigenvalue  7  is  complex,  with  the  imaginary  part  producing  the  modification 
of  the  streamwise  wavenumber  associated  with  the  change  in  frequency. 


t 


NUMERICAL  METHODS 


The  Falkner-Skan  equation  (2),  the  OSE  (9)  ,  and  the  secondary  dis¬ 
turbance  equations  (14)  are  numerically  solved. 

The  Falkner-Skan  profile  f(y)  is  calculated  using  a  4^  order  Runge- 
Kutta  integrator  coupled  with  a  shooting  scheme  to  satisfy  the  boundary 
condition  at  infinity.  An  accurate  value  of  /"(O),  necessary  to  initiate  the 
shooting  procedure,  is  obtained  from  a  parabolic  fit  of  f"(0)  as  function  of 
P  from  existing  data  [35],  thus  minimizing  the  number  of  iterations,  and, 
more  importantly,  avoiding  the  calculations  of  profiles  with  backflow  [35]  in 
adverse  pressure  conditions,  -0.198  <  P  <  0. 

For  the  fundamental  and  subharmonic  profiles  (18)  the  lowest  Fourier 
truncation  is  used:  v_2,v0,v2  for  the  fundamental,  and  v.^vj  for  the 
subharmonic.  The  dominance  of  the  first  terms  in  the  series  and  the  small 
error  of  the  truncation  have  been  proven  in  previous  works  [36],  [37],  Both 
the  OSE  and  the  secondary  disturbance  equations  are  solved  by  transform¬ 
ing  the  differential  equations  into  an  algebraic  form  through  the  use  of  a 
spectral  collocation  method  with  Chebyshev  polynomials  [38].  In  view  of 
the  finite  domain  of  definition  for  Chebyshev  polynomials,  the  semi-infinite 
domain  y  6  [0,co)  is  mapped  into  j?  G  [1,0)  by  the  transformation 

y  = - y0.  Letting  K  denote  the  number  of  collocation  points  used  in 

y 

the  analysis,  one  chooses  the  value  of  y0  such  that  half  the  collocation 
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points  yk  =  cos(  2A:— )  are  placed  inside  the  boundary  layer  displacement 
thickness,  wherein  the  disturbances  undergo  rapid  change. 

Without  loss  of  generality,  we  use  the  subharmonic  instability  problem 
to  discuss  in  detail  some  aspects  of  the  numerical  method.  We  approximate 
the  eigenfunctions  vn  =  (u„  ,  t/n)  in  series  expansion  (18)  by  a  finite  series 
of  odd  Chebyshev  polynomials, 

*»($)—  T,  ank^2k +  i(y)  >  ®»(y)  =  E  6^r2/  +  i(y)  (22) 

k  =  1  /  =  1 

The  use  of  only  odd  polynomials,  all  of  which  vanish  at  y  =  0,  causes  the 
boundary  conditions  at  that  end,  which  is  the  image  of  y  —  oo,  to  be 
automatically  satisfied.  Substituting  the  Chebyshev  series  representation  of 
the  velocity  components  into  the  governing  equations  and  solving  the  equa¬ 
tions  at  the  K  collocation  points  results  in  a  4 K  row  by  4K  +  6  column 
algebraic  system  for  the  unknown  coefficients  ank  and  b ni.  The  system  is 
made  solvable  (square  matrix)  by  the  addition  of  the  algebraic  form  of  the 
six  boundary  conditions  at  y  =  1,  which  is  the  image  of  the  wall  at  y  =  0, 


un  =  0  ,  vn  —  0  , 


For  elegance,  and  computational  simplicity,  the  ’’chain  rule  terms”  arising 
from  the  algebraic  mapping  are  incorporated  in  the  definition  of  the  deriva¬ 
tives  of  the  Chebyshev  polynomials.  The  ability  to  do  this  is  based,  clearly, 
on  the  linearity  of  the  equations;  for  example, 

dVn  ii“  dv~  _  1  _  j‘t  “  1  A 


dVn  y-  dV.  _  A 
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The  advantage  of  solving  the  problem  in  its  original  boundary  value 
format,  versus  solving  it  via  a  shooting  technique,  is  great.  A  global  eigen¬ 
value  search  procedure,  such  as  EISPACK  (Smith  et  al  [39])  can  be  applied 
to  the  algebraic  system  to  obtain  a  complete  spectrum  of  the  eigenvalues  at 
new  locations  in  the  multidimensional  parameter  space  (a,^,^,A).  Choos¬ 
ing  from  the  global  list  a  most  interesting  eigenvalue,  one  can  then  compute 
its  variation  during  a  local  walk  in  the  multidimensional  space.  A  Newton 
iteration  procedure  is  used,  which  computationally  less  demanding:  order 
(4K+6)3  versus  order  (4K+6)4  for  EISPACK. 

A  discussion  of  the  local  procedure  can  be  easily  done  by  introducing 
the  following  notation.  We  represent  the  algebraic  formulation  of  the  prob¬ 


lem  by 


G[  y]  =  0 


where  y  is  the  p-dimensional  vector  composed  of  the  unknown  coefficients, 


y —  (aik  >  ,  a. ik  >  &_i/) 


k  =  1  ,  K  +  1  ,  l  =  1  ,  K  +  2  ,  p  =  2(K  +  I)  +  2(K  +  2) 

and  G  :  Rp  — »RP  is  the  linear  operator,  acting  on  y,  constructed  by  the  spec¬ 
tral  collocation  method.  The  eigenvalue  a  appears  linearly  in  G,  while  7 
appears  up  to  fourth  power.  For  generality,  we  let  X  denote  either  of  the 
eigenvalues  and  indicate  the  functional  dependence  of  G  on  X  by  writing 


The  local  procedure  is  based  on  the  property  that  the  eigenvalue  X  and 
the  eigenvector  y  are  smooth  functions  of  the  parameters  (a  .  3  ,  R  ,  .4). 
Thus,  given  an  exact  pair  X0  ,  y0  at  (a 0  ,  3 0  ,  R0  ,  A0)  we  search  for 


X*  ,  y,  at  (ag  +  Aa  ,  0O  +  &p  ,  R0  +  &R  ,  A0  +  aA)  via  Newton  iteration. 
To  this  end,  we  augment  the  system  by  constructing 


x  =  (y  ,  X) 


x  €  Rp 


and  Fx  :  Rp  +  *— *•  Rp  +  1  is  obtained  by  adding  to  G  an  eigenvector  normali¬ 
zation  condition  (  for  example,  setting  the  last  coefficient,  b_n,  equal  to  1 
).  Note  that  F  is  not. a  linear  operator  when  X  represents  7,  i.e  in  a  spatial 
formulation.  Via  a  Taylor  series  expansion  we  write, 

0  =  F(xa)  +  v^x.[*,-x0]  +  0(  |x,-x0|2)  (28) 

where  VFX.  is  the  Jacobian  evaluated  at  x0,  acting  on  x,  -  xg.  Omitting 
the  higher  order  terms  we  rewrite  the  equation  in  the  form 

V^xr.,(xv-  XV-  1]  —  F{*v-  1)  (29) 

and  solve  for  xu  iteratively.  As  written,  the  Jacobian  is  updated  at  each 
iteration,  yielding  a  convergence  rate  of  second  order  at  a  computational 
cost  of  order  p3.  Suppressing  the  updating  to  selected  intervals  gives  a  con¬ 
vergence  rate  of  first  order,  but  decreases  the  computational  labor  to  order 
p2.  Since  20  to  30  collocation  points  are  necessary  for  good  accuracy,  the 
value  of  p  is  between  86  and  126,  so  that  the  selected  updating  of  the  Jaco¬ 
bian  can  greatly  reduce  the  numerical  work. 


RESULTS 


1.  Introduction 

In  low  turbulence  conditions,  such  as  those  that  exist  in  free  flight,  the 
subharmonic  type  of  secondary  disturbance  is  much  more  unstable  than  the 
fundamental  type  (Klebanoff  structure)  [40].  This  fact,  together  with  the 
need  of  extensive  computer  time  for  a  detailed  study,  led  us  to  focus  our 
attention  exclusively  on  the  subharmonic  modes.  The  study  of  pressure 
gradient  on  subharmonic  disturbances  is  complicated  by  the  functional 
dependence  of  the  eigenvalues  X  (growth  rates)  on  a  large  number  of 
parameters, 

X  =  X  (a  ,  0  ,  R  ,  A  ,  P)  (30) 

A  significant  portion  of  the  research  was  employed  for  orientation  in 
this  multi-dimensional  space.  Hence,  as  a  preparation  for  discussion  of  the 
results,  it  is  useful  to  mention  the  most  interesting  choices  of  parameter 
variation.  The  spanwise  wavenumber  selectivity,  given  by  the  variation  of 
growth  rate  with  0,  addresses  the  question  of  two  or  three-dimensionality  of 
the  disturbance.  This  variation  is  included  in  almost  all  the  graphs.  To  take 
into  account  the  change  of  the  nondimensionaiizing  variable  <5  with  stream- 
wise  location  we  replace  0  with  the  variable  b 

b  =  1000  0  /  R 

-  24  - 


i  >. 


which  describes  a  wave  of  fixed  physical  spanwise  wavenumber  as  it  travels 
downstream.  Other  important  parameters  are  the  TS  amplitude,  the  pres¬ 
sure  coefficient  P,  and  the  Reynolds  number.  The  influence  of  at  most  two 
parameters  can  be  clearly  represented  in  a  graph.  The  following  combina¬ 
tions  are  presented: 


ot  ,  R 

»  A 

fixed,  ; 

b 

,P 

varied. 

a  ,  R 

,P 

fixed  ; 

b 

,A 

varied. 

oc  ,  P 

,  A 

fixed  ; 

b 

,  R 

varied. 

0  ,  R 

,  A 

fixed  ; 

b 

,  oc 

varied. 

The  value  of  a  is  a  specified  parameter  in  temporal  analysis,  and  the  TS 
wave  frequency  u ;,  is  given  as  the  imaginary  part  of  the  eigenvalue  of  the 
OSE.  As  done  above,  we  define  a  new  frequency  ,  F,  independent  of  the 
downstream  change  of  the  nondime nsionalizing  variable  6, 

F  =  10®  w  /  R 

In  spatial  analysis,  the  value  of  the  frequency  F  is  specified,  and  the  TS 
wave  wavenumber  ar  is  given  as  the  real  part  of  the  eigenvalue  of  the  OSE. 
Hence  note  that  the  functional  dependence  on  o  in  (30)  is  changed  to  F 
when  dealing  with  results  from  the  spatial  analysis.  The  spatial  calculations 
were  conducted  mainly  at  F  =  124  and  F  —  83  to  match  experimental  con¬ 
ditions  [8]  ,[41] . 

2.  Effect  of  pressure  gradient 

We  begin  the  discussion  of  results  by  presenting  a  comparison  between 
theory  and  experiment  to  validate  the  accuracy  of  the  theoretical  model. 
Fig.  4  juxtaposes  experimental  [8]  and  theoretical  results  for  the  spatial 


amplitude  growth  of  TS  wave  and  subharmonic  mode  for  the  Blasius  profile. 
The  initial  TS  amplitude  is  =  0.44%  at  branch  I,  (R[  =  535)  and 
grows  to  1.23%  at  branch  II  (R[j  =  855).  A  virtual  leading-edge  correction 
of  20  cm  has  been  applied  to  the  experimental  data  in  order  to  bring  the 
experimental  Rir  =  840  into  agreement  with  the  theoretical  value,  855. 

The  initial  subharmonic  amplitude  was  chosen  at  B ^  =  0.00126%  in 
order  to  match  the  experimental  value  of  B  at  R //.  The  subharmonic 
amplitude  is  calculated  using 


R 


For  data  points  at  R  >  700  the  computed  growth  rate  closely  follows  the 
observed  slope,-  confirming  the  parametric  nature  of  subharmonic  instability 
in  the  initial  stages  of  growth.  The  measurements  at  R  <  700  were  taken 
shortly  downstream  of  the  vibrating  ribbon  {R  =575)  where  neither  the  TS 
wave  nor  the  subharmonic  mode  are  fully  established.  Fig.  5  shows  the 
selectivity  of  the  subharmonic  instability  mechamism  with  respect  to  6  at 
R  =700  and  853.  Based  on  the  maximum  spatial  growth  rate,  the  theory 
predicts  6  =  0.27  at  R  =  700  and  b  =  0.29  at  R  =  855.  However, 
natural  selection  of  b  in  experiments  also  depends  on  the  spectral  content  of 
the  background  noise.  The  observed  value  of  b  —  0.33  is  close  to  max¬ 
imum  amplification  throughout  the  range  of  measurement. 

Fig.  4  and  5  are  almost  indistinguishable  from  those  given  in  a  previous 
comparison  [2,  Fig.  9,  13  ]  in  which  the  growth  rates  are  computed  using  a 
temporal  formulation  and  converted  to  the  experimentally  observed  spatial 
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growth  rates  by  means  of  a  simple  transformation.  We  devote  an  entire  sec¬ 
tion  to  the  the  comparison  of  spatial  and  temporal  data. 

A  detailed  collection  of  experimental  data  for  transition  over  surfaces 
with  an  adverse  pressure  gradient  is,  today,  unavailable.  This  is  unfortunate 
because  the  anticipated  increase  in  the  error  due  to  the  shape  and  quasi¬ 
static  assumptions  cannot  be  assessed  by  comparison.  Indeed  there  is  evi¬ 
dence  [42]  that  near  separation  the  distortion  of  the  mean  flow  is  not  negli¬ 
gible.  However,  we  expect  the  current  results  to  capture  accurately  the 
quantitative  behavior  of  the  growth  rates. 

The  destabilizing  effect  of  an  adverse  pressure  gradient  on  the  dom¬ 
inant  subharmonic  mode  is  shown  in  Fig.  6.  A  decrease  in  the  pressure 
coefficient  increases  the  growth  rate  as  well  as  the  range  of  unstable  span- 
wise  wavenumbers.  A  likely  consequence  of  this  effect  is  the  amplification 
of  a  relatively  narrow  band  of  spanwise  wavenumbers  in  regions  of  favor¬ 
able  pressure  gradients,  which  sets  the  stage  for  the  downstream  develop¬ 
ment.  Most  surprising  in  Fig.  6  is  the  sharp  cutoff  at  b  ^  0.05  for  the 
inflectional  profiles,  indicating  strong  damping  of  the  vortex  pairing  mode  at 
a  TS  amplitude  of  1%.  Detuned  modes  (o  complex)  are  less  amplified  than 
the  dominant  mode  in  all  profiles  and  at  all  values  of  the  TS  amplitude  we 
have  investigated. 

While  the  adverse  pressure  gradient  enhances  the  growth  rates  of  both 
the  TS  wave  and  the  subharmonic  disturbance,  there  is  only  a  moderate 
effect  on  qualitative  features.  Dramatic  changes  occur,  however,  in  the 
eigenfunctions,  i.  e.  in  the  disturbance  velocity  profiles.  A  comparison  of 
mean  and  streamwise  rms  profiles  is  given  in  Fig.  7  for  the  Blasius  profile. 


P  =  0,  and  in  Fig.  8  for  the  Falkner-Skan  profile  near  separation, 


P  =  -0.18.  Fig.  7  shows  the  agreement  of  the  spatial  result  for  the  subhar¬ 
monic  mode  with  the  experimental  data  [8].  Comparison  with  [2,  Fig.  14] 
also  indicates  that  spatial  and  temporal  results  for  the  velocity  profiles  are 
hardly  distinguishable.  The  position  of  the  maxima  with  respect  to  the  criti¬ 
cal  layer  at  ye  (for  the  TS  wave  below,  and  for  the  subharmonic  mode 
above  ye)  is  qualitatively  unaffected  by  the  pressure  gradient.  However,  ye 
in  Fig.  8  is  further  away  from  the  wall.  A  second  maximum  develops  in  the 
TS  profile  owing  to  the  gradual  changeover  from  viscous  to  vortical  (inviscid 
inflectional)  instability.  Velocity  profiles  at  values  of  R  =  1168  and 
F  —  83  are  shown  in  Fig.  9.  We  note  by  comparison  to  the  profiles  at 
R  —  873  that  only  a  small  modification  is  brought  about  by  a  change  in  fre¬ 
quency  and  Reynolds  number.  Experimental  data  on  disturbance  velocities 
are  scarce.  At  the  present,  we  can  only  say  that  our  results  are  not  incon¬ 
sistent  with  measurements  [42,  Fig.  14]  in  the  separating  flow  over  an  air¬ 
foil. 

3.  Comparison  of  wall-bound  and  free  shear  layers 

Boundary  layer  profiles  near  the  point  of  separation  exhibit  high 
inflectionality  and  maximum  shear  stress  far  from  the  wall,  thus  bearing 
resemblance  to  a  free  shear  layer.  This  observation  leads  one  to  enquire 
under  which  conditions  the  two  types  of  flow  share  stability  characteristics. 
In  the  light  of  separation  and  reattachment  we  place  particular  attention  on 
the  preferred  spanwise  wavelength  of  subharmonic  modes,  in  other  words, 
on  the  two-dimensional  or  three-dimensional  nature  of  secondary  distur- 
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bances. 

The  inviscid  stability  of  a  spatially  periodic  free-shear  layer  has  been 
studied  numerically  by  Pierrehumbert  8c  Widnall  [20].  They  found  the  layer 
to  be  most  unstable  with  respect  to  a  two-dimensional  (0  =  0)  subharmonic 
mode,  and  the  growth  rate  monotonically  decreasing  with  increasing  0  [20, 
Fig.  5].  The  disturbances  undergo  pure  exponential  growth  and  are  tuned  (a 
real). 

Growth  rates  for  the  pairing  mode  and  a  three-dimensional  subhar¬ 
monic  mode  in  the  Falkner-Skan  profile  at  separation  are  shown  in  Fig.  10 
as  a  function  of  the  TS  amplitude.  Under  these  conditions,  three- 
dimensional  secondary  instability  is  dominant  up  to  TS  amplitudes  of  18%. 
At  low  amplitudes,  the  vorticity  concentrations  generated  by  the  TS  wave 
are  too  feeble  to  overcome  the  damping  effect  of  the  wall  on  the  pairing 
mode.  Vortex  pairing  takes  the  lead  at  very  large  amplitudes  of  the  primary 
disturbance,  but  still  competes  with  three-dimensional  phenomena.  This  is 
more  clearly  shown  by  the  growth  rates  as  a  function  of  the  spanwise 
wavenumber  b  in  Fig.  11.  There  is  a  strong  similarity  between  the  curve  for 
A  —  30%  and  the  curve  in  [20,  Fig.  5]  for  amplitudes  in  the  same  range. 
Note,  however,  that  the  neglection  of  nonlinear  terms  under  the  shape 
assumption  may  be  invalid  at  this  high  amplitude  level.  The  effect  of  the 
increase  of  vorticity  on  the  streamlines  can  be  seen  on  Fig.  12.  At  amplitude 
of  1%  the  vorticity  concentration  at  the  critical  layer  is  weak,  creating  a 
streamline  pattern  greatly  flattened  by  the  presence  of  the  wall.  At  ampli¬ 
tude  of  30%,  however,  the  effect  of  the  wall  is  not  pronounced  and  the 
streamlines  resemble  closely  those  of  a  free  vortex.  This  process  is  accom- 


panied  by  the  stability  changeover  to  a  basically  inviscid  vortical  mechanism. 

As  a  further  check  that  detuned  modes  are  always  less  amplified  than 
tuned  modes  a  comparison  of  growth  rates  was  made  for  different  values  of 
a.  For  the  separation  profile  (  P  —  -  0.198  ),  Fig.  13  shows  the  growth 
rate  of  the  most  unstable  fixed  mode  together  with  the  two  dimensional 
(0  =0)  running  mode  as  a  function  of  a.  For  the  most  unstable  fixed 
mode  0  varies  from  0.34  at  a  =  0.5,  to  0.08  at  a  =  0.05.  Clearly,  as  a  is 
decreased  from  the  upper  branch  to  the  lower  branch  of  the  TS  neutral 
curve  this  mode  always  dominates  in  growth  rate. 

Since  the  large  TS  amplitudes  required  for  the  pairing  mode  to  dom¬ 
inate  can  hardly  be  reached  by  viscous  growth,  we  can  conclude  that  within 
the  entire  attached  boundary  layer  the  development  of  secondary  instabili¬ 
ties,  and  hence  the  second  stage  of  transition,  will  be  dominated  by  three- 
dimensional  disturbances.  Through  this  result  it  is  fair  to  observe  that  the 
process  of  vortex  pairing  in  the  separated  shear  layer  is  controlled  not  by  the 
growth  of  vortex  pairing  prior  to  separation  but  by  the  amplitude  growth  of 
the  three-dimensional  disturbances.  This  observation  is  substantiated  by 
experiments  in  free  shear  layers  [43]  in  which  vortex  pairing  was  observed 
to  dominate  only  when  the  background  disturbances  where  sufficiently  low 
(hence,  low  amplitudes  of  the  three-dimensional  disturbances  in  the  boun¬ 
dary  layer  over  the  splitter  plate). 

4.  Temporal  to  spatial  transformation  for  secondary  disturbances 

In  view  of  the  reduced  computational  effort  needed  for  a  temporal 
eigenvalue  search,  as  opposed  to  a  spatial  search,  a  natural  interest  arises  in 


a  transformation  relating  the  two  types  of  growth  rates. 


Indeed,  one  may  even  hope  for  the  existence  of  a  simple  and  accurate 
relation  valid  over  a  wide  subdomain  of  the  parameter  space  since,  after  all, 
such  is  the  case  for  primary  (TS)  disturbances  (  Gaster’s  transformation 
[44].  For  a  transformation  in  non-parallel  flows,  see  [45]). 

The  transformation  for  secondary  disturbances  can  be  derived  in  a 
straight  forward  manner,  as  shown  below.  Unfortunately,  the  transforma¬ 
tion  takes  on  a  simple  form  only  in  regions  near  a  neutral  growth  point, 
where  both  growth  rates  are  small.  The  transformation  acquires  an  involved 
form  in  regions  of  large  growth  rates,  and  the  substantial  amount  of  compu¬ 
tation  required  may  undercut  its  usefulness. 

In  the  moving  coordinate  system  ( £,y,z,t )  the  velocity  field  solution  to 
the  secondary  disturbance  equations  takes  the  form 

g  V.  (»)«“<  (30) 

n  00 

where  a  and  7  are  the  complex  temporal  and  spatial  growth  rates,  respec¬ 
tively,  and  a  and  /?  are  the  streamwise  and  spanwise  wave  numbers,  respec¬ 
tively.  We  recall  that  solutions  with 

a  ~  a  t  +  ^il  represent  temporal  growth  m  both 
7  =  0  >  fixed  and  moving  coordinate  systems  (^^a) 


<7=7  c 
7  =  7  r  +  7 1 


represent  spatial  growth  in  fixed  coordinates 


(31b) 


The  disturbance  equations  functionally  relate  the  parameters  7, a. a,  and  3 
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F(~t,<r,oc,0)  =0  (32) 

provided  that  all  other  parameters  appearing  in  the  equations,  namely  R, 
frequency,  T-S  amplitude,  and  pressure  coefficient  P,  remain  constant. 

We  will  henceforth  assume  that  the  polysurface  described  by  equation 
F  is  smooth  and  the  partial  derivatives  with  respect  to  its  variables  exist  up 
to  the  order  required. 

Seeking  a  transformation  from  temporal  to  spatial  growth  we  re-write 
(32)  in  the  explicit  form 

T  =  /(r)  ,  r  =  (<7,a,0) 

We  choose  two  points  r4  and  rt  6  (ata,0),  and  connect  them  with  a 
smooth  path  c  (p)  parametrized  by  p,  with  end  points  c  (0)  =  rt  and 
c  (1)  —  t3  .  The  values  of  7  at  the  two  end  points  of  the  path  are  related 
by  the  first  fundamental  theorem  of  calculus, 

1 

7(0  -  l{rt)  =  f  grad  (  F)  •  (33) 

This  is  a  most  general  transformation.  In  the  special  case  when  ra  and 
r(  correspond  to  spatial  and  temporal  growth  rates  respectively,  the  end 
points  are  according  to  (31), 

ra  =  (7  °c,a,0)  and  rt  ==  {a°,ot,0) 

where  7°  is  the  spatial  and  a0  the  temporal  growth  rate.  The  path  c(p) 
now  keeps  the  value  of  a  and  0  constant,  and  (33)  reduces  to 


l’c 


where  use  has  been  made  of  the  Cauchy  -  Riemann  condition 
dq  _  d~ir  .  <hi  _  £21  _  •  £2 1 

9(7  9(7  r  9(7  r  <9(7  ,•  9(7  ,• 


Equation  (34)  is  the  desired  transformation.  Note  that  the  spatial  eigen¬ 
value,  which  is  the  unknown,  appears  on  both  sides  of  the  transformation. 
Next  we  want  to  reduce  equation  (34)  to  a  form  involving  only  data 
obtained  at,  or  near,  the  endpoint  c  (0)  corresponding  to  temporal  growth. 
Therefore,  we  expand  the  integrand  in  a  Taylor  series  about  <7  ? 

-  ».  +  ».'(»,  -  of)  +  tfi*'  -  o?)S  +  -  •  (36) 


where 


etc.  (37) 


Inserting  the  Taylor  series  into  (5)  we  obtain  the  following  approximate 
transformations  for  real  eigenvalues, 


One  term  :  7  ®  =  J*  g0  da  r 


(38a) 


(c-  J-) 
90 


Two  term  :  7®  =  J  |?0  +  -  ^°')}  r 


where 


a2  i;  +  aiir  +  fl5  =0 


a2  =  ?'ac2 

«i  =  [9 0  ~  -  1 


(38b) 


■Slf 

... 


i 


A  three-term  Taylor  series  yields  a  cubic  algebraic  equation  for  7®,  a  four- 
term  series  yields  a  quartic,  and  so  on.  Since  approxim  .'.ons  of  higher 
order  require  numerous  values  of  7  =  l{&)  near,  and  at,  a  °  for  the  calcu¬ 
lation  of  higher  derivatives  and  also  require  numerical  solution  to  the  associ¬ 
ated  algebraic  equation,  high  order  approximations  become  too  cumbersome 
to  calculate.  The  results  for  a  one-term  and  trwo-term  transformation  are 
given  below.  Blasius  flow  Re  =  1168  ,  F  —  83  ,  a  =  0.276  ,  A  '=="  1% 
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Clearly,  for  0  near  the  neutral  stability  point,  both  the  one-term  and 
two-term  approximation  yield  good  results.  However,  away  from  this  point 
the  approximation  becomes  inaccurate.  The  error  due  to  truncation  is  large 
at  higher  values  of  7®  and  <7®,  where  the  function  7  —  ~i{cr)  exhibits 
strong  curvature.  This  can  be  seen  in  Fig.  14.  The  function  7  —  7(a)  is 
quasi  linear  at  0  =  0.09  while  for  higher  values  of  0  the  curvature  under¬ 
goes  strong  bending  as  it  approaches  the  spatial  endpoint  7  ®c. 

Under  strong  curvature  numerous  terms  in  the  Taylor  series  expansion 
(34)  must  be  kept  for  accuracy,  and  due  to  the  accompanying  increase  of 
computational  work  the  transform  loses  its  labor  saving  advantage.  Expan¬ 
sions  other  than  a  Taylor  series  could  yield  better  convergence. 

The  plot  of  7  =  ~t(0)  at  fixed  Re,ar,F,  and  T-S  amplitude  is  shown  in 
Fig.  15,  along  with  the  transformed  values  at  three  0  locations.  The  zeroth 
order  approximation 


is  the  one  used  in  previous  works  [2],  and  gives  the  error  upper  bound  for 
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CONCLUSION 


The  preseat  Investigation  finds  that,  under  free-flight  conditions  (low- 
turbulence)  ,  the  secondary  disturbances  in  boundary  layers  with  favorable 
and  adverse  pressure  gradients  (up  to  separation)  behave  similarly  to  those 
in  the  Blasius  boundary  layer,  with  the  most  noticeable  effect  of  the  pres¬ 
sure  gradient  being  a  decrease  (favorable)  or  increase  (adverse)  of  the 
disturbance’s  growth  rate.  In  particular,  the  dominant  disturbances  are 
always  three-dimensional  and  tuned  (real  eigenvalues)  to  the  phase  speed  of 
the  TS  wave.  Only  at  experimentally  unreachable  amplitudes  in  the  neigh¬ 
borhood  of  20-30%  does  the  secondary  stability  of  the  boundary  layer 
become  similar  to  that  of  a  free  shear  layer,  whereat  the  vortex  pairing 
dominates  in  growth.  This  fact  is  of  importance,  since  it  is  conjectured  that 
the  rapid  evolution  of  three-dimensional  disturbances  prior  to  the  separation 
of  the  boundary  layer  may  cause  the  separated  shear  layer  to  break  down 
into  wall-bound  turbulence. 

A  comparison  of  growth  rates  obtained  from  a  temporal  analysis  and  a 
spatial  analysis  shows  that  the  results  from  the  two  formulations  can  be 
related  easily,  with  an  acceptable  amount  of  error,  but  an  involved  transfor¬ 
mation  is  necessary  for  an  accurate  change.  Thus  temporal  analysis  remains 
a  useful  tool  due  to  its  lower  computational  demand,  and  spatial  analysis  is 
of  use  when  an  accurate  comparison  with  experimental  data  is  needed. 
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Due  to  a  lack  of  experimental  data  for  transition  in  boundary  layers 
with  non-constant  pressure  gradient,  a  much  needed  comparison  of  theory 
and  experiment  is  not  possible.  Good  agreement  in  the  Blasius  boundary 
layer,  however,  suggests  that  at  least  qualitatively  the  presented  results  are 


APPENDIX 


Herein  we  substitute  the  TS  streamfunction  and  the  normal  mode 
form  of  the  secondary  disturbance  (19)  into  the  governing  equations 
(15a,b),  and  arrive  at  the  final  formulation  of  the  problem.  First,  we  divide 
the  differential  equations  into  a  sum  of  operators  acting  separately  on  u  and 


Equation  (15a), 

T  /  \  [l  o  trr  \  d  d  )  „  ,  d2U0  dv 

L<e>“  HP  ‘  {u-'  c)a7-  Ttrv  +  —a? 

+  2.8** 

J  dx^dx2  dy2  dz2>  dxdy  dxdy 

,  n  dy2-?/;  du  ,  d2y2^ 

+  2~aT  +  “~i?- 

P-(*„)  _  (*LL.  S±±.)^v  +  2 £±J£l.  +  j!jLi£z 

'  dx  dy  dy  dx'  <9X2  dxdy  dxdy''  dy2 

_  d2v  +  d2v  >  dy2,0  dv  dy2-^  du  ^d2^2j) 

dx 2  dz 2  dx  dy  dy  dx  dxdy 


Equation  (15b), 


M*  ( n )  = 


Mu  ( v )  = 


Q*  ( 0» » )  = 


f  1  O  f  rr  sd  d  |,  d2U  d2U  V 

IT7'  -  (c,“-  c)aT-  aTl(i^  +  ^-> 

1  2  _  /■  rr  N_d_  d  d2v  jg,  d2v 

V  0  dx  dt  ,  dxdy  dy  dz 2 

,  d’d)  d d^  d d2t/>  w  d2u  d2a .  d2T/>  d2u 

dx  dy  dy  dx  dxdy  Qx 2  dr2  ^  dx2  dxdy 


q-(«,»)  =  i*LL .  i±l - <£L.)J£± l 

V  dx  dy  dy  dx  dxdy  dxdy 


d2d>  d2v  _  d2U>  d 
dx2  dy 2  dy2  d 
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The  general  form  (no  truncation)  for  the  TS  wave  streamfunction  is 

=  A  £  »,<■*'  (40) 

m  =  -  oo 

The  general  form  (no  truncation)  for  the  secondary  disturbance  equation  is 

v3(x,y,z,t)  =  c,0*cffte”,x  £  vn(y)emdz  (41) 

»  =  -  00 

where  a  =  a/2  and  v„  =  (  un  ,  vn).  Introducing  (40)  and  (41)  into  equa¬ 
tions  (15a,b)  and  collecting  terms  multiplied  by  exponentials  with  like  index 
yields  the  following  infinite  sets  of  equations, 

Equation  (15a), 

L(«„)  t!4  j  [P,(e„,«J)  +  !**(*,..#/)]  =  0  (42a) 

m  =*-  oo 

Equation  (15b), 

oo 

M*(«n)  +  M’(«„)  +  RA  £  [Q“(0m,«;-)  +  <Ww,u;)]=0  (42b) 

m  =»-  oo 

where  j  —  n  -  2m 

The  infinite  order  system  (42a,b)  must  be  truncated  for  computational  rea¬ 
sons.  Under  the  shape  assumption,  we  truncate  the  TS  streamfunction  (40) 
to  two  terms;  4>_l  and  <t>lf  and  for  the  secondary  disturbance  the  lowest 
truncation  is  used;  v_2,va,v2  for  the  K-type  disturbance  (fundamental), 
and  v-i>vi  for  the  subharmonic  disturbance.  The  system  (42a,b)  separates 
into  two  uncoupled  equations  for  the  K-type  and  the  subharmonic  distur¬ 
bance: 


■r. 

-J  - 
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EQUATIONS  FOR  K-TYPE 


Equation  (42a), 

L(u2)  +  flAP“(^,us  )  +  RAPv(<j>l  ,v0  )  =  0 
L(t;0)  +  RAPu{4>vu_2)  +  RAPB{<i>l,v_2). 

+  RA  P*(^_i,«2)  +  RAPv(<t>_  vv2)  =  0 
L(t/_2)  +  RAP'(<j>_vu0)  +  RAPv(<t>_vv0)  —  0 

Equation  (42b), 

Ma(«2)  +  Mv(  v2)  +  RA  Q“(<^i  ,u0  )  +  RAQHi^V'  )  =  0 
M“K)  +  Mv(v0)  +  RAQr{<f>lfu_2)  +  RAQv{<j>l,v_2). 

+  RA  Qm(<t>_lfu2)  +  RAQ?{4_ltv2)  =  0 
M“(«_2)  +  Mw(u_2)  +  RAQu(<i>_l,u0)  +  RAQv(<p_vva)  —  0 


EQUATIONS  FOR  SUBHARMONIC 


Equation  (42a), 

L(ui)  +  RAPu(4>l,u_l)  +  RAPv(<t>l,v_1)  —  0 
L(v.,)  +  fiAP^^^j)  +  RAF'{<i>_l,vl)  =  0 

Equation  (42b), 


M“(  «i)  +  Mv(  ut)  +  RAQm^u.J  +  RAQ'it^v.i)  =  0 
M“(«_!)  +  M*(  v_  t)  +  RA  Q“(0_  ! ,  u  t )  +  5AQ*(<^_  pyj  =  0 
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Figure  1.  Effect  of  two-dimensional  and  three- 
dimensional  disturbances  on  a  separating  shear  layer. 
Top,  vortex  pairing  and  permanent  separation.  Bottom, 
three-dimensional  vortex  stretching,  and  breakdown  into 
wall-bound  turbulence. 
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Figure  3.  Construction  of  the  periodic  basic  flow.  The 
weak  streamwise  variation  of  boundary-layer  thickness 
and  TS  wave  amplitude  are  neglected.  The  local 
conditions  at  X0  are  extended  into  a  periodic  flow. 
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Figure  5.  Spatial  growth  rate  7  versus  spanwise 
wavenumber  b  for  the  conditions  of  Fig.  4.  The  TS 
amplitudes  are  A  =  0.74%  at  R  =  700,  .4  =  1.23%  at 
R  —  855.  The  vertical  line  marks  the  observed 
wavenumber  b  =  0.33  . 


Figure  10.  Temporal  growth  rate  a  versus  TS  amplitude 
A  for  tuned  modes  with  6  =  0  (a)  and  6  =  0.36  (b). 
P  —  -  .198,  R  =  1000,  a  =  0.5,  F  =  245. 


Figure  11.  Temporal  growth  rate  a  versus  spanwise 
wavenumber  6  for  various  TS  amplitudes  A  in  a  profile 
at  separation.  P  =  -  .198,  R  —  1000,  a  =  0.5, 
F  =  245. 


Figure  13.  Most  unstable  tuned  mode  and  detuned 
vortex  pairing  as  function  of  a.  P  =-  -0.198  R  =  1000 
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(ABSTRACT) 

The  transition  from  laminar  to  turbulent  flow  in  boundary- layers  occurs 
in  three  stages:  onset  of  two-dimensional  TS  waves,  onset  of  three- 
dimensional  secondary  disturbances  of  fundamental  or  subharmonic  type, 
and  onset  of  the  turbulent  regime.  In  free  flight  conditions,  subharmonic 
disturbances  are  the  most  amplified. 

Recent  modeling  of  the  subharmonic  disturbance  as  a  parametric  insta¬ 
bility  arising  from  the  presence  of  a  finite  amplitude  TS  wave  has  given 
results  in  quantitative  agreement  with  experiments  conducted  in  a  Blasius 
boundary- layer.  The  present  work  extends  the  analysis  to  the  Falkner-Skan 
family  of  profiles,  and  develops  a  formulation  for  spatially  growing 
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disturbances  to  exactly  match  the  experimental  observations. 


Results  show  that  subharmonic  disturbances  in  Falkner-Skan  flows 
behave  similarly  to  those  in  a  Blasius  flow.  The  most  noticeable  effect  of  the 
pressure  gradient  is  a  decrease  (favorable)  or  an  increase  (adverse)  of  the 
disturbance’s  growth  rate.  Due  to  the  lack  of  experimental  data,  a  com¬ 
parison  of  subharmonic  growth  rates  from  theory  and  experiment  is  limited 
to  the  Blasius  boundary-layer.  A  comparison  of  results  from  the  spatial  for¬ 
mulation  with  those  previously  obtained  from  a  temporal  formulation  shows 
the  difference  to  be  small.  A  connection  between  disturbance  growth  in  a 
separating  boundary-layer  profile  and  a  free  shear  layer  is  presented.  A 
modification  of  Gaster’s  transformation  from  temporal  to  spatial  growth 
rates  for  secondary  disturbances  is  given. 


